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Abstract 

The thermal correction to the energy of Casimir-Polder interaction of atoms with a suspended 
graphene membrane described by the Dirac model is investigated. We show that a major impact 
on the thermal correction is made by the size of the gap in the energy spectrum of graphene quasi- 
particles. Specifically, if the temperature is much smaller than the gap parameter (alternatively, 
larger or of the order of the gap parameter), the thermal correction is shown to be relatively small 
(alternatively, large). We have calculated the free energy of the thermal Casimir-Polder interac- 
tion of atoms of He*, Na, Rb, and Cs with graphene described by both the hydrodynamic and 
Dirac models. It is shown that in exact computations using the Dirac model, one should use the 
polarization operator at nonzero temperature. The computational results for the Casimir-Polder 
free energy obtained in the framework of hydrodynamic model of graphene are several times larger 
than in the Dirac model within the separation region below 2 ^m. We conclude that the theoretical 
predictions following from the two models can be reliably discriminated in experiments on quantum 
reflection of different atoms on graphene. 

PACS numbers: 31.30.jh, 34.35. +a, 12.20.-ni, 42.50. Ct 



I. INTRODUCTION 



The Casimir-Polder force acting between a microparticle and a macrobody is caused by 
the existence of quantum fluctuations of the electromagnetic field. This phenomenon is a 
manifestation of the more general dispersion forces at relatively large separations where the 
retardation of the electromagnetic interaction becomes significant (at separations below a 
few nanometers the electromagnetic interaction can be considered as instantaneous and the 
dispersion force is usually referred to as the van der Waals force). Quantum theories of the 
van der Waals and Casimir-Polder forces were developed by London [l| and (for the case of 
an atom interacting with an ideal metal plate) by Casimir and Polder ^ , respectively. For 
an atom interacting with a plate made of some real material described by the frequency- 
dependent dielectric permittivity, the unified theory of the van der Waals and Casimir-Polder 
forces was developed by Lifshitz js] (see also monographs 4-6|). 

Nowadays, the Casimir-Polder interaction between different atoms and material plates 
(cavity walls) has attracted much attention in experiments on quantum reflection [t-Q]. 
This is a process in which an ultracold atom is reflected by an attractive atom-wall poten- 
tial or, in other words, it is an above-barrier reflection. It was shown [10] that quantum 
reflection is very sensitive to the specific form of the Casimir-Polder interaction. Theoret- 
ically, calculations of the Casimir-Polder forces between various atoms and plates made of 



different materials were performed 



11 



14| on the basis of the Lifshitz theory. The results ob- 



;ained were compared [15| with those calculated using a simple phenomenological potential 
8|,Q,0. 

Recently, special attention has been directed to the Casimir-Polder interaction with car- 
bon nanostructures, such as graphene, carbon nanotubes, and fuUerenes 18|. Computations 



were performed using the phenomenological density-functional approach 19N22|. second- 



order perturbation theory 



23| and, for multiwalled carbon nanotubes with at least several 



walls, using the Lifshitz theory [2j] . For one-atom-thick carbon nanostructures, the concept 
of dielectric permittivity is not immediately applicable (it becomes applicable, for exam- 



ple, to thin fullerene fi' 



ms and can be used to deduce the dynamic polarizabilities of single 



fullerene molecules [25|). In order to extend the Lifshitz theory to this case, the reflection 



coefficients o: 
model 1261 



electromagnetic oscillations on graphene were found using the hydrodynamic 



2i\ . For this purpose, graphene was treated as an infinitesimally thin positively 
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charged flat sheet, earring a homogeneous fluid with some mass and negative charge densi- 
ties. The obtained reflection coefficients depend on temperature only through the Matsubara 
frequencies. Calculations of the Casimir- Polder interaction in the framework of the hydro- 

3 (see also ^^). 



dynamic model are presented in Refs. 



28 



The hydrodynamic model of graphene is only a crude approximation. It does not take into 
consideration that low-energy excitations in graphene are massless Dirac fermions except for 



the fact that they move with a Fermi velocity, rather than with the speed of light [18|, |31|. 
From this, it follows that at low energies the dispersion relation for quasiparticles is linear 
with respect to the momentum measured relative to the corner of the graphene Brillouin 
zone. The Dirac model of graphene assumes that the dispersion relation is linear at any en- 
ergy. Using this assumption, the corresponding reflection coefficients for the electromagnetic 
oscillations at zero temperature have been found 32| different from those obtained with the 
hydrodynamic model. The obtained coefficients depend on the polarization operator in an 
external electromagnetic field calculated in the one-loop approximation in three-dimensional 
spacetime. Note that the exact value of the gap parameter A of quasiparticle excitations 
entering the polarization operator is not known. In Ref. 33j it was shown that Dirac model 
leads to much smaller values of the van der Waals (Casimir-Polder) atom- graphene inter- 
action than the hydrodynamic model at separations below 100 nm. Keepirig in mind the 



precision achieved in experiments on quantum refiection of ultracold atoms [7|-l9|, [34, |35| , it 
was concluded that it is possible to experimentally distinguish between the predictions of 
the Dirac and hydrodynamic models of graphene. 

The thermal effect in the Casimir and Casimir-Polder interactions is a subject of debate 
[g], [3^. For the thermal interaction between two graphene sheets it was argued 36 1 that 
the thermal regime starts at rather low temperatures T (or, respectively, short separation 
distances of tens of nanometers at room temperature) because the value of the effective 
temperature is determined by the Fermi velocity vp rather than the speed of light c. This 
conclusion was, strictly speaking, obtained in a nonretarded regime. A fully consistent 
quantum-field version of the Dirac model at nonzero temperature was presented in Ref. js^] ■ 
In that paper the temperature-dependent polarization operator and respective coefficients 
for electromagnetic oscillations on graphene at nonzero temperature were derived. Compu- 
tations were performed for a graphene layer with the gap parameter equal to zero interacting 
with a metallic plate. For this system the conclusion of Ref. 36| on the existence of large 
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thermal corrections at short separation distances was quahtatively confirmed (with a proviso 
that the characteristic distance separating zero- and high-temperature regimes depends on 



the fine-structure constant rather than on the Fermi velocity). In Ref. [38| quantum refiec- 
tion of ultracold atoms from thin films, semiconductor heterostructures and graphene was 
considered at T 7^ with graphene being described by both the hydrodynamic and Dirac 
models. It was concluded that suspended graphene membranes produce higher quantum- 



refiection probabilities than bulk matter. The remark made earlier [33|] that already achieved 
experimental precision allows discrimination between the predictions of different models of 
graphene was supported. This conclusion was obtained using the polarization operator at 
zero temperature 32|. Keeping in mind that suspended graphene membranes with rather 



large diameter (of 55 /im) are already available 39|, it seems pertinent to perform a full 
quantum-field theoretical investigation of the thermal Casimir-Polder force between various 
atoms and such membranes in the framework of different models of graphene proposed in 
the literature. 

In this paper, we calculate the free energy of the thermal Casimir-Polder interaction be- 
tween atoms of metastable helium (He*), sodium (Na), rubidium (Rb), and cesium (Cs) and 
a suspended graphene sheet. Graphene is described using either the hydrodynamic model 
or the Dirac model and the obtained results are compared. In the framework of the Dirac 
model, the polarization operator at nonzero temperature is employed and the dependences 
on the value of the gap parameter A are investigated. We demonstrate that under the con- 
dition ksT ^ A, where is the Boltzmann constant, the thermal correction is relatively 
small, but if A < ksT there are large thermal corrections to the Casimir-Polder force. At the 
same time, the magnitude of the thermal correction strongly depends on the atom-graphene 
separation. The value of the Fermi velocity is shown to be of no crucial influence on the 
magnitude of thermal correction for atom-graphene system. The novelty of this paper is 
that the interaction of different atoms with graphene was investigated using the full Dirac 
model at nonzero temperature with the temperature-dependent polarization operator. In so 
doing it was found that the relative size of thermal correction strongly depends on the gap 
parameter. Our computations show that experiments on quantum reflection are capable to 
discriminate between the predictions of the hydrodynamic and Dirac models of graphene. 
One can also conclude that in the framework of the Dirac model it is necessary to take 
into account the dependence of the polarization operator on temperature. According to our 



4 



results, the use of the zero-temperature polarization operator in computations at T = 300 K 
can lead to large errors in theoretical predictions for some values of parameters. 

The paper is organized as follows. In Sec. II we present a brief formulation of the two 
models of graphene, introduce the reflection coefficients and illustrate the limiting case of 
zero temperature for the Dirac model. Section III contains the investigation of dependences 
of the thermal correction to the Casimir-Polder energy on the gap parameter. In Sec. IV 
the distance dependences of the Casimir-Polder interaction of different atoms used in exper- 
iments on quantum reflection with graphene are calculated. In Sec. V the reader will find 
our conclusions and discussion. 



II. TWO DIFFERENT MODELS FOR THE REFLECTION COEFFICIENTS ON 
GRAPHENE 

The unified expression for the van der Waals and Casimir-Polder free energy of an atom 



Ml. 



This formula can 



interacting with a planar structure is given by the Lifshitz formula 
be expressed in terms of refiection coefficients of the electromagnetic oscillations on this 
structure (in our case on graphene) in the following way: 



X {2yVTM««, y) - Cf Km(«0, y) + rTE««, y)]} ■ 

Here, a is the separation distance between the atom and the graphene sheet, a{iC,i) is the 
dynamic polarizability of an atom calculated along the imaginary Matsubara frequencies = 
2TikBTl/h with / = 0, 1, 2, . . . , the dimensionless Matsubara frequencies are Q = (.i/oJc, the 
characteristic frequency is defined as Uc = c/(2a), and ttm.te are the refiection coefficients 
for two independent polarizations of the electromagnetic field (transverse magnetic and 
transverse electric, respectively). The prime near the summation sign multiplies the term 
with Z = by 1/2. Note that the dimensionless frequencies Q are functions of the separation 
a. 



5 



A. Hydrodynamic and Dirac models of graphene 



As was discussed in Sec. I, there exist two different models of reflection coefficients for 
graphene. In the framework of the hydrodynamic model the reflection coefficients take the 
form 



26 



Ky + Cf 

rTE {Ki ,y) = {Kuy) = - -J— ■ (2) 

K + y 

Here, K = 2aK and K = 6.75 x 10^ m^^ is the characteristic wave number of the graphene 
sheet which corresponds to the frequency uk = cK = 2.02 x lO^^rad/s. The dimensionless 
variable y is connected with the projection of the wave vector on the graphene sheet k± by 
the equation y = 2a{k\ + /c^Y^"^. As was mentioned in Sec. I, the hydrodynamic model 
does not take into consideration that the dispersion relation for quasiparticles in graphene 
is linear with respect to momentum. It should be mentioned also that the parameter K in 
the hydrodynamic reflection coefficients ([2]) is temperature- independent. 

In the framework of the Dirac model the reflection coefficients are expressed in terms of 
the components of the polarization tensor js^] 



r-TM«;, y) = -~ — -— — (3) 
{y^ - C^)fltr - y^tloo 



(y2_c2)(n,, + 2y)-y2n, 



00 



where the dimensionless components Iloo.ir are connected with the dimensional Hoo^tr by the 
equation noo,tr(^0;Z/) = (2a/^)noo,tr(^0) !/)• The explicit expressions for the components of 
the polarization operator at nonzero temperature were obtained in Ref. js^]. In terms of our 
dimensionless variables the component IIoo can be identically represented in the following 
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form: 

noo««,y) 



- Ci ) / dx 



x{l — x) 



A2 + x(l-x)/(0,z/) 



n 1/2 + ^2 / dx 



Jo 



(4) 



A2 + Cfx(l - x) 



cos(27r/x) + e 



A2 + x(l-x)/(0,z/) 



cosh g{T,(i,y) + cos(27r/x) 



Here, a = /{he) ~ 1/137 is the fine-structure constant, A = A/{JtjjJc) is the dimensionless 
gap parameter, r = ATTakBT/{hc) = Q/l, vp = vp/c is the dimensionless Fermi velocity, the 
chemical potential fi is assumed to be equal to zero, and the dimensionless functions / and 
g are defined as 



fiCi,y) = vly' + il-dl)Cf 



(5) 



27V 



A^ + xii-x)fiCi,y) 



1/2 



In a similar way, for the sum of two spatial components of the polarization tensor 11^^ in our 
dimensionless variables one obtains 



Utri^Cl,y)=8a[y^ + f{Cl, 



dx- 



x{\ — x) 



A2 + x(l-x)/(0,y) 



1/2 



+ - 



JO 



dx (6) 



X < — In [1 + 2 cos(27r/x)e-5("'^"^) + e'^^^^'^'^y^ 



;i-2x)- 



sinf27rZx) 



2 cosh g{T, (i, y) + cos(27r/a;) 

A2 + x(l - x)[(l - ^2)2^2 _ -4^2] cos(27r/s) + e-^^^^^^'V^ 



, /I \\^^'^ cosh g(r, 6, y) + cos(27rZx) 

A2 + x(l -x)/(0,y)J '^"^^ ^ ' 

It is easily seen that in the limiting case of zero temperature (T — )■ 0, r — )■ 0) we have 
g ^ oo and the components of the polarization operator (jlj) and become much simpler 



noo«,Z/) = a 
^tr{iC,y) = a 



fiCy) 

y' + fiCv);, 



^oo{<,y), 
^oo{i(,y), 



(7) 



where 

x(l — x) 



$oo«,2/) = 8v//(C,l/) / dx- -^'^ (8) 



and C is now the continuous dimensionless frequency. Calculating the integral in Eq. ([8]) we 
arrive at 

<loo«,y) = 4A + 2v//(C,l/) 



A2 

1-4- 



arctanVSS^l. (9) 
2A ^ ^ 



7(C,l/)_ 

After the substitution of Eq. ([7]) in Eq. ([3]) the reflection coefficients at zero temperature 
take the form 

a?/l>(C,y) 



ay<l>(C,y) + 2/(C,2/) 



a<l>(C,y) + 22/ 



These equations were obtained in Ref. [32 1. 



B. Properties of reflection coefficients in the Dirac model 

Now we return to the consideration of general reflection coefficients ([3]) in the Dirac model 
with the polarization operator and (|6]) defined at any nonzero temperature. As can be 
seen from Eq. ([3]), for y = Q we have 

rTM{iChy)\y^^^ = I, rTE«/,Z/)|j^^^^ = -1, (11) 

as it holds for any Q and y for ideal metal plane at both zero and nonzero temperature. At 
zero temperature, for a graphene sheet described by the Dirac model, Eq. ( II ip . however, 
does not hold because in accordance with Eq. ([71) the quantity Iloo{i(,y)\^^^ = and both 
reflection coefficients in ([3]) become indeterminate form (i.e. zero/zero). The reflection 
coefficients at T = are converted to a determinate form in Eq. (fTOj) . It is seen that the 
coefficients ttm.te f|TO|) at y = ( are not equal to the limiting values of the coefficients (ITT]) 
when T — )■ (i.e., not equal to 1 and -1, respectively). This means that in the Dirac model 
the reflection coefficients calculated under the condition y = d are discontinuous functions 
of temperature at the point T = (note that unlike the hydrodynamic model, the reflection 
coefficients of the Dirac model depend on T not only through the Matsubara frequencies 
but also explicitly through the polarization operator). 
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To illustrate the behavior of the reflection coeflicients in the Dirac model, Fig. 1(a) shows 
by the solid lines rTM,TE(^Ci) u) as functions of y > (i = 0.49 at T = 300 K, a = 300 nm and 
A = (rxM is positive and tte is negative). As can be seen in Fig. 1(a), ttm decreases with 
decreasing y and abruptly jumps to unity in the vicinity of y = (i. In Fig. 1(b) the same 
is shown in an enlarged scale. The dashed line in Fig. 1(b) shows the coefficient rTM(^Ci)l/) 
calculated with the polarization operator at zero temperature [i.e., by substituting IloQ{i(i,y) 
from Eq. ([7]) in Eq. (|3]) instead of using the operator at T = 300 K defined in Eq. (jl])]. From 
Fig. 1(b) it is seen that at / = 1 the use of the polarization operator ta T = leads to almost 
the same values of the TM refiection coefficient as the use of the operator at T = 300 K (the 
relative difference between the solid and dashed lines for almost all values of y is of about 
4%). Significant difference between the two calculation methods arises only within a very 
narrow interval from (i to (i + 10"''. For Z > 2 all differences under discussion become even 
smaller. 

The solid line illustrating the behavior of the refiection coefficient rTE(^Ci5 Z/) ^s a function 
of y in Fig. 1(a) almost coincides with the horizontal coordinate axis in the scales used. In 
the close proximity of y = (i, the refiection coefficient abruptly jumps to minus unity. On 
an enlarged scale the behavior of rTE(^Ci)2/) as a function of y is shown in Fig. 1(c). From 
the comparison of Fig. 1(b) and Fig. 1(c), one can conclude that \rTE{i(i,y\ ^ |'"tm(^Ci)Z/| 
for all y with exception of only a very narrow vicinity of the point y = Qi- Note also that the 
values of rTE(^Ci;l/) calculated using the polarization operators at zero temperature and at 
T = 300 K are indistinguishable in the scales of both Fig. 1(a) and 1(c). The same holds for 
I > 2. Thus, the use of the polarization operator ((Tj) instead of (jlj) and in calculations 
of rTE(^Ci; ?ziwith / > 1 leads to even smaller errors than for r^j^uiKi^y)- 



In Ref. 



37| it was proposed to use the refiection coefficients ffTOl) taken at zero temper- 



ature in all terms of the Lifshitz formula with / > 1 and restrict the application of the 
exact refiection coefficients ([3]), (jlj) and (j6j) to only the zero- frequency term / = 0. In our 
computations performed below we determine the accuracy of this prescription. 
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III. DEPENDENCE OF THE THERMAL CORRECTION ON A GAP PARAME- 
TER 



In this section we calculate the free energy of thermal Casimir-Polder atom-graphene 
interaction and the thermal correction to the Casimir-Polder energy in the framework of the 
Dirac model of graphene with different values of the gap parameter A. The exact value of 
A is yet unknown. The upper bound on A is approximately equal to 0.1 eV, but it might be 



also much smaller 



18| . As an atomic system interacting with graphene, we choose an atom 



of metastable helium He* often used in experiments on quantum reflection 17| . To perform 
computations using the Lifshitz formula ([T]), one needs the dynamic atomic polarizability of 
He* as a function of the imaginary frequency. For many atoms the dynamic polarizability 
can be expressed with sufficient precision, using the single-oscillator model 

,. a(0) 
"^'"'^'^ = l + (u;gK)Cr ^''^ 
where a(0) is the static polarizability and cjq is the characteristic absorption frequency. 

Specifically, for He* we have a{0) = a^'^ (0) = 315.63 a. u. = 4.678 x 10~^^m^ (where one 
atomic unit of polarizability is equal to 1.482 x 10~^^ m^) and uq = u^"^ = 1.18 eV = 
1.794 X lO^^rad/s 4^. Note that the use of highly accurate dynamic atomic polarizabilities 



(see, for instance, the polarizability of He* determined with a relative error 10 



4lj) lead 



to only small relative deviations from the results obtained using Eq. (fT2l) . For example, 
for He* atom near an Au wall these deviations decrease from 3.9% at a = 10 nm to 0.03% 
at a = 1 fim jlll . fisl ]. Our computations show that for graphene the contribution of the 
term with / = in Eq. ([T]) is dominant even at short separations, i.e., the contribution of 
the static atomic polarizability is of most importance. Because of this, for graphene the 
single-oscillator model leads to even more accurate results than for metallic walls. 



A. Casimir-Polder free energy as a function of temperature 

We have substituted Eqs. (l3])-(l6l) and ( fT2l) in Eq. ([T]) and performed computations of 
the Casimir-Polder free energy J-" as a function of temperature at atom-graphene separation 
a = 1 fim for the values of a gap parameter equal to A = 0.1, 0.05, 0.01 eV and for A = 0. 
The computational results for the quantity are presented in Fig. 2, where the four 

lines from the lowest to highest correspond to the decreasing values of A (the lowest line 
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is for A = 0.1 eV). As can be seen in Fig. 2, the characteristic behavior of the free energy 
differs significantly for different gap parameters. At T = the values of the Casimir-Polder 
energy E{a) = J^{a,0) depend heavily on A, whereas at T = 300 K there is only a slight 
dependence of the Casimir-Polder free energy J^{a,T) on A. This allows estimation of A 
from the comparison between the measurement data and computational results in the region 
of moderate temperatures from 100 to 150 K. Furthermore, the larger is A, the wider is the 
temperature region where J-" remains constant with the increase of temperature. In such 
temperature regions the thermal correction to the Casimir-Polder energy is negligibly small. 
Below we discuss this point in more detail. 

To check quantitatively an accuracy of the prescription [s^l that in all terms of the 
Lifshitz formula with / > 1 one can use the polarization operator at T = 0, we repeated 
the same computation as above, but this time with the operators @ and ([6]) for / = and 
([7]) for all / > 1. The obtained results cannot be distinguished visually from those shown 
in Fig. 2. The largest deviations between the computational results obtained using different 
calculation procedures (with the polarization operators found at T 7^ for all / and for only 
/ = 0) hold with the gap parameter A = 0. In this case, the magnitudes of the free energy 
obtained using the prescription are smaller than those in full computations by 0.6%, 0.3%, 
0.1%, and 0.06% at T = 300 K and separations a = 50, 200, 500, and 1000 nm, respectively. 
At a = 3 um the relative difference between the computational results obtained using the 

n 

two procedures is as small as 0.002%. One can conclude that the prescription of Ref. [37| 
leads to very accurate results and can be used in subsequent computations. 



B. Thermal correction as a function of separation 

Now we calculate the thermal correction to the Casimir-Polder energy in the interaction 
of He* atom with graphene as a function of separation. The relative thermal correction at 
a temperature T is defined as 

The computations are performed at T = 300 K using Eq. ([1]) with full polarization operator 
for / = and zero-temperature operators for / > 1. The computational results within the 
separation region from 10 to 500 nm are shown in Fig. 3(a) by the six solid lines from the 
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highest to the lowest corresponding to the values of the gap parameter A = 0.1, 0.05, 0.025, 
0.01, 0.001, and OeV, respectively. Note that the two lowest lines obtained for two smallest 
values of the gap parameter are almost coinciding. It is interesting that at T = 300 K the 
magnitude of the relative thermal correction at each separation increases monotonically with 
increasing A. The largest thermal correction is achieved for A = 0.1 eV. At T = 300 K the 
same properties hold at shorter separation distances below 100 nm. This separation region 
is ahown in an enlarged scale in Fig. 3(b). 

As suggested by Fig. 2, the monotonous dependence of the relative thermal correction on 
A is not universal and does not hold at any temperature. As an example we have computed 
the thermal correction to the energy of He*-graphene interaction at T = 100 K. This is the 
maximum temperature until which the free energy is nearly constant for A = 0.1 eV. The 
computational results within the separation region from 10 to 500 nm are shown in Fig. 4(a) 
by the three solid lines from the lowest to the highest corresponding to the following values 
of the gap parameter: A = 0.1, 0.05, and 0.025 eV. As is seen in Fig. 4(a), for A = 0.1 eV the 
thermal correction is very small over the entire separation region. This is in accordance with 
the computational results shown in Fig. 2. The computational results for the same thermal 
correction at T = 100 K, but with the values of the gap parameter equal to A = 0.025, 0.01, 
0.001, and OeV, are shown in Fig. 4(b) by the solid lines from the highest to the lowest line, 
respectively. For illustration purposes, the line with A = 0.025 eV is reproduced in both 
Fig. 4(a) and 4(b). From the comparison of Fig. 4(a) with Fig. 4(b) one can see that at 
T = 100 K the dependence of the thermal correction on A is not monotonous. From Fig. 4(b) 
it can be seen also that there is a noticeable difference between the thermal corrections for a 
graphene with A = 0.001 eV and with A = 0. Similar results are obtained for other atoms. 

One can conclude that the size of thermal correction to the Casimir-Polder interaction of 
an atom with a graphene sheet depends essentially on the size of the gap in the spectrum of 
graphene quasiparticles. From our computations it follows that if the inequahty ksT <^ A 
is satisfied with a large safety margin, then the relative thermal correction is small (note 
that T = 100 K corresponds to ksT = 0.0083 eV to be compared with A = 0.1 eV). On the 
contrary, if A < ksT, then the thermal correction is large (in so doing, the thermal correction 
can be also large for ksT smaller but not much smaller than A). Keeping in mind that the 
value of A for graphene is not yet known [isl, the predictions of the Dirac model concerning 
the size of thermal correction in atom-graphene interaction remain uncertain. 
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IV. DISTANCE DEPENDENCE FOR THE CASIMIR-POLDER INTERACTION 
OF DIFFERENT ATOMS WITH GRAPHENE 

We come now to the Casimir- Polder interaction of graphene with atoms of He*, Na, 
Rb, and Cs at room temperature T = 300 K but at different separation distances, as is of 
interest for experiments on quantum reflection. Here, we compare the computational results 
obtained using both the hydrodynamic and Dirac models for graphene. 

A. Atom of metastable helium 

We start with an atom of He* and compute the Casimir-Polder free energy ([1]) with the 
reflection coefficients ([2]) of the hydrodynamic model. The computational results for a^\J^\ 
as a function of separation from 50 nm to 5 /im are presented in Fig. 5(a) by the dashed line. 
Then we repeated computations using Eq. ([T]), but with the Dirac reflection coefficients Q, 
which contain the polarization operator (E]) obtained at T = 300 K [as was noted in 
Sec. II, it is sufficient to use this operator only with / = term and replace it by a more 
simple operator ([7]) in all terms with / > 1]. In this case the computational results for the 
gap parameter in the region from to 0.01 eV are shown by the solid line. Note that even 
the use of larger gap parameters up to 0.1 eV leads to almost imperceptible shift of the solid 
line in Fig. 5(a). This is explained by the fact that at T = 300 K there are only minor 
differences between the Casimir-Polder free energies computed with different A (see Fig. 2). 

In Fig. 5 (a) we also plot by the dotted line the computational results obtained using 
the Lifshitz formula ([1]) and the Dirac model with A = at zero temperature [i.e., using 
the polarization operator ([7]) with all Matsubara terms with / > 0]. As can be seen from 
Fig. 5(a), the computational results shown by the solid line are by factors of 1.27, 1.29, and 
1.29 larger than those shown by the dotted line at separations 1, 3, and 5 ;um, respectively. 
This confirms that in order to perform precise computations at T = 300 K, one should use the 
polarization operator at the same temperature (at least in the zero-frequency contribution of 
the Lifshitz formula). From Fig. 5(a) it is also seen that the predictions of the hydrodynamic 
model at the separations 0.2, 0.5, 1.0, 1.5, 2, and 3 yum (the dashed line) are by factors of 
3.83, 2.96, 2.07, 1.62, 1.36, and 1.12 larger than the predictions of the Dirac model (the solid 
line). This allows reliable discrimination between theoretical predictions of the two models 
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by the measurement data of experiments on quantum reflection. (Note that in Ref. |33| . 
where computations were performed at a < 100 nm using the polarization operator at T = 0, 
the predictions of the Dirac model were overestimated by approximately a factor of 1.5 due 
to an error in the computer program. As a result, the differences between the predictions of 
two models at short separations were underestimated.) 

From the comparison of solid and dotted lines in Fig. 5(a), we have already found errors 
arising from the use of polarization operator at T = in all terms of the Lifshitz formula 
(an underestimation of \J-'\ by approximately a factor of 1.3). This conclusion was obtained, 
however, from the zero-temperature polarization operator with A = (as was demonstrated 
above, computational results with the T-dependent polarization operator are not sensitive to 
the value of A at T = 300 K). To investigate the role of the gap parameter in computations 
using the polarization operator at T = 0, in Fig. 5(b) we plot J-"! as a function of separation 
computed with the full operator (jlj), (the solid line) and with the operator ([7]) for 
A = 0, 0.001, and 0.01 eV (the dotted lines from the highest to the lowest, respectively). 
Note that the solid and the highest dotted lines reproduce the respective lines in Fig. 5(a) 
at separations from 50 nm to 1 /im. As can be seen from Fig. 5(b), for nonzero A the 
underestimation of \J^\ when using the polarization operator at T = is much larger than 
for A = 0. Thus, at a = 1 fim this underestimation is by the factors of 1.83 and 6.2 for 
the gap parameter A = 0.001 and 0.01 eV, respectively. This adds importance to the use 
of full temperature-dependent polarization operator (jl]), ([6]) in computations performed for 
subsequent comparison with the experimental data. 

B. Atoms of sodium, rubidium and cesium 

The Casimir- Polder free energy for other atoms used in experiments on quantum reflection 
can be computed in a similar way. For Na the dynamic polarizability can be presented 
by Eq. ([T2Dwith a^^(O) = 162.68 a.u. = 2.411 x lO-^Sm^ and co^^ = 2.14 eV = 3.25 x 
10^^ rad/s 42|. The computational results for the quantity a"^! J-"! for Na-graphene interaction 
at T = 300 K are presented in Fig. 6(a) as function of separation by the dashed line (the 
hydrodynamic model), the solid line [the Dirac model with the T-dependent polarization 
operator (j4]), ([6])], and by the dotted line [the Dirac model with the polarization operator 
(I7j) at T = 0, A = 0]. It can be seen that qualitatively the computational results in Fig. 6(a) 
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are similar to those in Fig. 5(a) for a He* atom. In the case of Na, however, the values 
of a'^\J^\ are smaller than for He* at all respective separations. The difference between 
the predictions of the hydrodynamic and Dirac models can be easily discriminated by the 
measurement data of experiments on quantum reflection. Thus, at separations 0.2, 0.5, 1.0, 
1.5, and 2 /im the predictions of the hydrodynamic model are larger than the predictions 
of the Dirac model with the T-dependent polarization operator by factors of 4.04, 3.05, 
2.10, 1.63, and 1.37, respectively. The use of the polarization operator at T = again 
underestimates the predictions of the Dirac model (by the factors of 1.27, 1.29, and 1.29 at 
separations a = 1, 3, and 5 /xm, respectively). 

In Figs. 6(b) and 6(c) similar results are presented for atoms Rb and Cs interacting with 
graphene. To perform computations of the Casimir-Polder free energy using different models 
of graphene, we have used the single-oscillator model f|T2l) with the following parameters: 
aRb(^Q) ^ 319.9 a.u. = 4.73 x lO-^^m^, = 5.46 eV = 8.3 x 10^^ rad/sksl and a^'^(O) = 
403.6 a.u. = 5.981 x lO-^^m^ and u^^ = 1.55 eV = 2.36 x lO^^rad/s 44|. As can be 
seen from the comparison of Fig. 6(b) with Fig. 6(a) and Fig. 6(c) with Fig. 6(b), the 
magnitudes of the free energy for Rb are larger than for Na and the respective magnitudes 
for Cs are larger than for Rb at all separation distances. This is explained by the fact that 
the static polarizability of Cs is the largest one. From Figs. 6(b) and 6(c) one can observe 
the same pattern between the dashed, solid and dotted lines, obtained for the hydrodynamic 
model and for the Dirac model with T-dependent and T-independent polarization operators, 
respectively, as was discussed on the basis of Figs. 5(a) and 6(a). 

In Figs. 5(a) and 6(a-c), the distinction between the predictions of the hydrodynamic 
model (the dashed lines) and the Dirac model with T-dependent polarization operator (the 



solid lines) deserves special 
same qualitative behavior 



attention. The point is that the dashed lines demonstrate the 



151] as the Casimir-Polder free energy of atoms interacting with 



walls made of ordinary real materials, such as Au or Si. By contrast, the Dirac model for 
graphene predicts quite a different behavior of the free energy, where it is nearly inverse 
proportional to the third power of separation. Keeping in mind that at separations below 
2 fim the predictions of both models may differ by a factor of 2 and even up to a factor 
of 4, it seems appropriate to raise a question which model is in better agreement with 
the experimental data. This question can be answered through experiments on quantum 
reflection of different atoms on suspended graphene membranes. 
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V. CONCLUSIONS AND DISCUSSION 



In this work we have investigated the Casimir-Polder interaction of different atoms with 
graphene in the thermal regime. Special attention was paid to the Dirac model of graphene 
and to the conditions when thermal effects become essential. We have confirmed the con- 
clusion js^] that sufficient precision can be achieved by using the temperature-dependent 
polarization operator in the zero-frequency term of the Lifshitz formula alone and calculat- 
ing all the other terms using the polarization operator at T = 0. We have found that a 
major impact on the thermal correction (i.e. whether it is small or large) is made by the 
size of the gap A in the spectrum of graphene nanoparticles. According to our results, if 
the condition fc^T ^ A is satisfied with a large safety margin, the thermal correction to 
the Casimir-Polder interaction of an atom with graphene is relatively small. In future it is 
planed to obtain the asymptotic behavior of the free energy under this condition analytically 
and to verify the fulfilment of the Nernst heat theorem in the Dirac model. If, however, 
A < ksT, the thermal correction is relatively large. 

Furthermore, we have calculated the free energy of the Casimir-Polder interaction of dif- 
ferent atoms used in experiments on quantum reflection (such as He*, Na, Rb, and Cs) with 
a suspended graphene membrane. All computations were performed for room temperature 
within a wide separation region in the framework of two models for graphene proposed in 
the literature: the hydrodynamic and Dirac ones. It is important to note that both these 
models are approximations and a priori it is not possible to decide which of them provides 
a better theoretical description of the Casimir-Polder force (for example, the Drude model 
of metals describes correctly the electric properties of metals and their dielectric response 
to external electromagnetic fields, but is in disagreement with the experimental data on 
measurements of the Casimir force between metallic test bodies jo, Isol). 




Our computational results allow to conclude that the hydrodynamic and Dirac models of 
graphene lead to both qualitatively and quantitatively different results for the free energy 
of atom-graphene Casimir-Polder interaction. On the qualitative side, we have arrived at 
quite different dependences of the free energy on separation distance. Quantitatively, the 
magnitudes of the free energies computed using the two models differ by a factor of 2 
and even by a factor of 4 at different separations below 2 fim. A difference in theoretical 
predictions being as large as by this amount assures a reliable discrimination between the 
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hydro dynamic and Dirac models of graphene by comparing with the experimental data 
on quantum reflection. Keeping in mind that suspended graphene membranes have been 



already produced 
future. 



391 ]. one could expect the resolution of this problem in the immediate 
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FIG. 1: (Color online) (a) The reflection coefficients ttm.te on graphene calculated at ^ = Ci, 
T = 300 K and a = 300 nm as functions of the dimensionless variable y using the polarization 
tensor at T = 300 K are shown by the positive- and negative- valued solid lines, respectively, (b) 
The solid and dashed lines show ttm computed using the polarization tensor at T = 300 K and 
T = 0, respectively, (c) The solid line shows tte- The scales in (b) and (c) are enlarged as 
compared to (a). 
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FIG. 2: (Color online) The Casimir-Polder free energy of He*-graphene interaction at a = 1 //m 
multiplied by the fourth power of separation is shown as a function of temperature by the four lines 
from the lowest to the highest for the gap parameter A = 0.1, 0.05, 0.01, and OeV, respectively. 
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FIG. 3: (a) The relative thermal correction to the Casimir-Polder energy of He*-graphene interac- 
tion at T = 300 K multiplied by the fourth power of separation is shown as a function of separation 
by the solid lines from the highest to the lowest for the gap parameter A = 0.1, 0.05, 0.025, 0.01, 
0.001, and OeV, respectively, (b) The same is shown at separations from 10 to 100 nm. 
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FIG. 4: The relative thermal correction to the Casimir-Polder energy of He*-graphene interaction 
at r = 100 K multiplied by the fourth power of separation is shown as a function of separation by 
the solid lines (a) from the lowest to the highest for the gap parameter A = 0.1, 0.05, and 0.025 eV 
and (b) from the highest to the lowest for the gap parameter A = 0.025, 0.01, 0.001, and OeV. 
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FIG. 5: (Color online) The Casimir-Polder free energy of He*-graphene interaction at T = 300 K 
multiplied by the fourth power of separation is shown as a function of separation (a) by the 
dashed, solid and dotted lines using the hydrodynamic model, the Dirac model with T-dependent 
polarization operator, and with the polarization operator at zero temperature, respectively, (b) 
The same quantity is shown where the dotted lines from the highest to the lowest are computed 
using the polarization operator at zero temperature with the gap parameter A = 0, 0.001, and 
0.01 eV, respectively. 
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FIG. 6: (Color online) The Casimir-Polder free energy of atom-graphene interaction at T = 300 K 
multiplied by the fourth power of separation is shown as a function of separation by the dashed, solid 
and dotted lines using the hydrodynamic model, the Dirac model with T-dependent polarization 
operator, and with the polarization operator at zero temperature, respectively, for atoms of (a) 
Na, (b) Rb, and (c) Cs. 
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